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Abstract 

Recent progress in hadron structure calculations within lattice QCD is reviewed. Results on key 
observables such as the axial charge, the quark momentum fraction and the spin content of the 
nucleon are discussed with focus on open issues. Lattice QCD studies of the j*N ^ A transition 
as well as the A form factors are also presented. 



1 Introduction 

The recent progress in the numerical simulation of the fundamental theory of the strong interactions, 
Quantum Chromodynamics (QCD), has been remarkable. Improvement in algorithms coupled with 
increase in computational power have enabled simulations to be carried out at near physical parameters 
of the theory. This opens up exciting possibilities for an ab initio calculation of experimentally measured 
quantities as well as for predicting quantities that are not easily accessible to experiment. During the 
last decade, results from simulations of QCD have emerged that already provide essential input for a 
wide range of strong interaction phenomena as, for example: i) The QCD phase diagram [U [21 El S! 
relevant for quark-gluon plasma searches at RHIC and LHC that probe the structure of our universe at 
t ~ 10~^^s after the big-bang at temperatures T ~ 10^^ C; ii) The structure of hadrons [51I61I71I81 flOlITT] 
that formed in our universe at t ~ 10~^ s. Key hadronic properties such as the nucleon axial charge, 
the quark momentum distribution and spin content of the nucleon as well as form factors of other 
hadrons, resonances and decays are being investigated within the experimental programs of major 
facilities such as JLab, LHC and Mainz, iii) Nuclear forces [13] that determine the large scale structure 
of the universe t ~ 10^ years after its birth. Although lattice QCD calculations of three- and four- 
baryon systems appeared recently, it is estimated that exa-scale computing resources will be required 
to extract accurate results and to extend to studies of light nuclei within lattice QCD |14j . 

In this work we will focus on hadron structure calculations [12] using state-of-the art lattice QCD 
simulations. After presenting some results that highlight the progress made in the meson sector, we will 
discuss the evaluation of key observables that probe the structure of the nucleon such as the form factors 
(FFs) and moments of generalized parton distributions (GPDs). Understanding nucleon structure from 

* Lecture presented at the Erice School on Nuclear Physics: From Quarks and Gluons to Hadrons and Nuclei, organized 
by A. Faessler and J. Wambach. 
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first principles is considered a milestones of liadronic pliysics and measurements of tfie electromagnetic 
nucleon form factors were first carried out more tlian 50 years ago. Despite tlieir long history of 
measurements, recent double polarization experiments greatly improved the accuracy and revealed 
unexpected features, namely the ratio of the proton electric to magnetic form factor, jJ'pG'^^^q^) j G\j{q^) , 
instead of being approximately constant, falls off almost linearly with the momentum transfer squared, 

[15]. This behavior is conjectured to be due to two-photon exchange terms that were previously 
neglected from the analysis [16J and it motivated new dedicated experiments to measure these form 
factors to higher precision and larger momentum transfers [I^. Compared to the electromagnetic 
(EM) form factors, the nucleon FFs connected to the axial- vector current are more difficult to measure 
and therefore less accurately known. An exception is the nucleon axial charge, qa-, which is precisely 
determined from neutron /3-decay and provides a benchmark quantity for lattice QCD techniques. 

Consequently we will present lattice QCD results on the 7*A^ — >■ A transition. The 7*A^ — )■ A is 
well studied experimentally as a probe of nucleon/ A deformation. In contrast to the nucleon and the 
7*A^ — > A FFs, the EM FFs of the A have not been measured, except from the A magnetic moment. 
Therefore they provide an ideal example of observables where lattice QCD can make predictions. We 
review these calculations as well as their implication on the shape of the A. Our study of the A and 

to A systems include the corresponding matrix elements of the axial-vector current that provide 
information on the axial couplings used in chiral expansions. 



2 Introduction to Lattice QCD 



Before we explain how the relevant hadronic matrix elements are determined we briefly outline the 
lattice formalism that enables their extraction. 



4j(n)4 



The starting point is a definition of the theory on a four- 
dimensional Euclidean space-time lattice. The lattice acts as a 
non-perturbative regularization scheme with the lattice spac- 
ing a providing an ultraviolet cutoff at highest allowed mo- 
mentum vr/a. Gauge fields are defined as links U^{n) between 
adjacent lattice sites and quarks are defined at each lattice 
site as anticommuting Grassmann variables belonging to the 

fundamental representation of SU(3). 
U^(n)=eig3AM(n) 

Using these fundamental quark and gluons degrees of freedom one constructs an appropriate action 
such that when a — )■ (and the lattice volume— >■ oo) one recovers the continuum theory. The con- 
struction of the appropriate operators with their renormalization is then carried out in order to extract 
physical quantities. This discrete formulation of QCD in Euclidean time is referred to as lattice QCD 
(LQCD). It can be simulated on the computer using methods analogous to those used in Statistical Me- 
chanics allowing calculation of matrix elements of any operator between hadronic states in terms of the 
fundamental quark and gluon degrees of freedom. LQCD, therefore, provides a well-defined approach 
to calculate observables non-perturbatively starting directly from the QCD Langragian with only input 
parameters the coupling constant and the quark masses. 

In what follows we will show results obtained on an isotropic hypercubic grid taking the spatial 
distance between lattice sites as to be the same as the temporal one a-r i.e. a = as = ax- We will also 
limit ourselves to zero temperature taking the spatial size Ls less than the temporal size Lt- Anisotropic 
lattices where as > at are also being used mainly for the study of excited states [HI [19]. Although 
LQCD provides an ab initio calculation of hadronic properties, the discretization of space-time and the 
numerical simulation on a finite volume introduce artifacts that may lead to systematic errors, which 
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must be carefully investigated before comparing to experimental results. Of particular relevance, for 
the observables discussed in this work, are the following issues: i) Finite Volume: One needs to perform 
simulations on different volumes to study finite volume effects. Volume studies of observables considered 
here, have shown that finite effects are small for lattice sizes such that Lsm.^ ~ 3.5. A consequence of 
having a finite volume is that only discrete values of momentum are allowed. For periodic boundary 
conditions (b.c.) the momenta allowed are in units of 2'k/Ls- Twisted b.c. have been used to obtain 
observables at small non-zero momenta, ii) Finite lattice spacing: LQCD simulations for at least three 
values of the lattice spacing are needed in order to extrapolate results to the continuum limit, iii) 
Renormalization constants: Computation of the appropriate renormalization constants is needed in 
order to relate lattice matrix elements to continuum results, iv) Heavier than physical pion masses: 
Current LQCD results are obtained using simulations that are typically performed at heavier than 
physical pion masses requiring chiral extrapolation. However, simulations with pions at their physical 
mass have been achieved [20] and more are foreseen in the near future, v) Fourier transforms: One 
needs to Fourier transform lattice results computed in coordinate space numerically. For large values 
of momentum transfer, results become noisy and typically only values of momentum transfer squared 
= — ~ 2 GeV^ can be achieved. 



2.1 Lattice techniques 

In this section we briefly explain how hadron masses and matrix elements of local operators are com- 
puted. Like the continuum theory, the QCD lattice action can be written a.s S = Sg + Sp, where Sg 
contains a purely gluonic part, written in terms of the gauge links [/^, and Sp contains the kinetic 
energy of the quarks and the interaction terms. The fermionic action is bi-linear in the quark fields 
ip^x) i.e. Sp = J2n,j '4'{T^)Dnji^{j), where n and j denote the lattice sites. The exact form of the lattice 
Dirac matrix, D, depends on the discretization scheme used for the fermions. Within the Feynman 
path integral representation the vacuum expectation value of a gauge invariant operator B is given by 

^ / d[u]d[m^] B[U, V;, ^ij]e-''.m-sm^,i^] 

After integrating over the fermionic degrees of freedom, one obtains 

< n\B\n >=\ j d[U] det(L'[[/])B[[/,D-i[[/]]e-^«[^l , Z = j d[U] det{D[U])e-^o^^^ , (2) 

in which the fermionic determinant appears and we obtain a factor of DJ^\U] for each possible Wick 
contraction of fermion pairs —ipnipj- One performs the path integrals numerically by stochastically 
generating a representative ensemble of gauge configurations {U} according to the probability 
exp {— + In {det{D[U]))} /Z given in Eq. ^ and then averaging over these gauge configurations: 

1 ^ 

<fiisifi>= hm . (3) 

The time consuming part of a LQCD calculation related to hadron properties is the generation of an 
ensemble of gauge configurations and the computation of the inverse of the fermionic matrix D, which 
yields the quark propagatoiQ In many applications only a column of is required. However, in 
calculating e.g. the isoscalar nucleon FFs the full (spatial) inverse is needed as these involve diagrams 
where an external probe couples to a sea quark. Calculation of these fermionic loops is therefore much 
more difficult and results on the observables presented in this work have not taken into account these 
contributions. Historically, LQCD simulations were done in the quenched approximation that sets 
det{D) = 1 in Eq. ^ leaving a local action Sg[U] making such simulations much easier. Nowadays, 
however, unquenched simulations are performed, which include the det{D). 



^In multi-baryon systems the contractions needed for the computation of observables is also a time-consuming step [H] , 
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Figure 1: Simulation cost of TMF using Ls = 
2.1 fm, a = 0.089 fm function of the pion 
mass to the p-meson mass [21j. The physical 
point is showed by the dotted horizontal line. 



The cost of these simulations can be estimate by 
using the cost formula: 
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where the coefficients c^, cl and depend on the 
type of discretization used for the fermionic ac- 
tion. State-of-the-art simulations use improved al- 
gorithms that take advantage of the mass precon- 
ditioner trick ^22j and use multiple time scales in 
the molecular dynamics updates. These improve- 
ments reduce the required simulation time making 
simulations at the physical value of the pion mass 
(physical point) feasible. The cost for generating 
10^ independent Nf = 2 twisted mass fermions 
(TMF) gauge configurations as the pion mass de- 
creases at fixed a and Ls is shown in Fig. [TJ One 
finds Cm ~ 4 in contrast to ~ 6 that charac- 
terized Wilson fermion simulations ten years ago. 
Based on current simulations the cost at the phys- 
ical point is estimated to be 0{1) Teraflop-year. 



2.2 Recent results 

2.2.1 Spectrum of low-lying baryons 

Masses of low-lying hadrons are extracted from the vacuum expectation value of two-point functions: 

G(t,p1 = Ee-^''^^(^|rJ„Tx^(f,t)x^(o,o) \n), (4) 

X 

with the projection matrix = (1 -|- 7^)/2. The interpolating fields, Xh{x), are operators in the 
Heisenberg representation that create a trial state with the quantum numbers of the hadron h that we 
want to study. The spectral decomposition of the two-point function, in the infinite lattice size limit. 
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Figure 2: Nucleon effective mass using local and 
smeared interpolating fields. 



Figure 3: Ratio of nucleon three-point to two-point 
functions for two sink-source separations func- 
tion of the current insertion time {t — ti)/a. 



4 



can be written as 



G{t,p) = Y.m\Xk\nm'e 



—mot 



(5) 



where En{p} is the energy of the n^^ state. For single particle states En{p) = ym'^+fp. By \h > 
we denote the lowest eigenstate of QCD with the quantum numbers of Xh with mass mo obtained by 
setting p = in Eq. ^ and overlap 2^(0) = | (i7|x/i|/i(0)) p with the trial state For large time 

separation t between the source and the sink the unknown overlap factor |(fi|x/i|/i(0))p and exponential 
time dependence cancel in the ratio mee{t) = — log G{t,0)/G{t — 1,0) , which therefore becomes time 
independent (plateau region) and can be fitted to a constant to yield mo- This is demonstrated in Fig. [2| 
where we show the nucleon effective mass as a function oit/a obtained using local interpolating fields 
and smeared ones. Smearing creates a quark field with support at several lattice sites instead of at one. 
It is a technique used to optimize the overlap Zh with the ground state achieving faster ground state 
dominance as shown in Fig. [2j For the numerical evaluation of such two-point functions one needs quark 
propagators from the fixed source (0, 0) to all spatial x-sink points for several time separations. This 
requires the computation of only one column for each spin/color component of the quark propagator, 
namely f (x, t) = D^^{x, t; 0, 0) obtained by solving the equation D{y, x)v{x) = 6{y — (0, 0)). 
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Figure 4: The low-lying baryon spectrum com- 
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TMF 



2 + 1 Clover [23j and Nf = 2 



In Fig. |4] we show recent results on the low-lying 
baryon spectrum obtained using two types of Wil- 
son improved fermion action Sp- The BMW Collab- 
oration uses a Clover term in Sp and smeared gauge 
links. Simulations were performed with two degen- 
erate u- and d-quarks and a strange quark fixed 
to its physical mass {Nf = 2 + 1) using the mass 
of the Q. The ETM Collaboration uses a twisted 
mass term in Sp, which provides automatic 0{a) 
improvement, and Nf = 2. The strange valence 
quark mass was fixed using the kaon mass. Both 
collaborations analyzed configurations generated at 
3 lattice spacings: a = 0.125,0.085,0.065 fm de- 
termined by the S-mass in the case of BMW and 
a = 0.089, 0.070, 0.056 fm, set by the nucleon mass 
in the case of ETMC and extrapolated the results to 
the continuum limit and to the physical pion mass. 
One observes that the results using different dis- 
cretization schemes are in agreement and that both 
reproduce the experimental values. This is a signif- 
icant validation of LQCD techniques. 



2.2.2 Hadron form factors 

Calculation of hadron matrix elements is more involved and requires the evaluation of a three-point 
function, depicted schematically in Figs. |5] and |6| The corresponding expression for the nucleon three- 
point function is 

G,,{r,q,t) = e'^'-HQ\r^^x:{xf,tf)O,,{x,t)x^^{0)\Q) . (6) 

Xf,X 

We first comment on the evaluation of the connected diagram. From Eq. (|6| it is clear that two spatial 
sums are involved, one over the spatial coordinates of the operator and one over the spatial coordinates of 
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the final state. This means that the propagator from the operator insertion, denoted by x in Figs. [5] and 
|6| to the sink, denoted by Xf, involves all spatial columns of D^^. The trick to evaluate this inverse is to 
perform the sum over by solving the equation a;)V(a;) = S{xf) — > V(x) = J2xf D~^{x,Xf)S{xf) 

with an appropriately constructed source S{xf) that combines the two quark propagators from the fixed 
source at with the hadron state at tf. The solution V{x) is the sequential (backward) propagator from 
the sink to the operator with the summation over Xf done automatically. Using the symmetries of the 
Dirac operator the forward sequential propagator can be easily constructed. This so called 'fixed sink 
method' is used in most recent calculations of three-point functions and takes its name from the fact 
that the quantum numbers of the hadron state at the sink enter into the construction of the sequential 
propagator and must therefore be fixed. It also requires fixing the sink-source time separation tf — ti, 
final momentum and spin projection matrix F. Inserting the operator, which can be done at all values 
of X, and summing over with the appropriate Fourier phase and propagator starting at ti and ending at 
t yields the connected three-point function. Therefore, within this scheme and for each non-degenerate 
quark flavor, two inversions involving the Dirac matrix enable one to compute the three-point function 
for all possible momentum transfers q and operators O. 

To extract the nucleon matrix element {N{pf)\0\N{pi)) we study the large Euclidean time behavior 
of an appropriately defined ratio of the three-point function and two-point functions ^12j given by 



G{pi,tf-t)G{0,t)G{0,t 



f) t»i 



G{0,tf) \G{0,tf-t)G{pi,t)G{pi,t 



n^^(g,F) , (7) 



f) 



where we have set pf = 0. Like in the case of the effective mass, this ratio is defined so that the unknown 
overlaps Zn{p) as well as the time dependence arising from the time evolution cancel, yielding a time 
independent quantity (plateau), which signals identification of the nucleon state from the tower of QCD 
states with the same quantum numbers as the nucleon. Fitting to this plateau value we can extract the 
matrix element {N [p f)\0\N [pi)) and from this, depending on the choice of (9, the FFs or moments of 
GPDs. However, the identification of the plateau region is much more delicate as compared to the case 
of the effective mass. This is demonstrated in Fig. |3] where we show results for two sink-source time 
separations. As can be seen, while one might think that there is already a plateau for the smaller time 
separation, increasing the sink-source separation changes the value of the plateau, which in turn means 
that excited states still contribute significantly. For the nucleon form factors and pion masses larger 
than about 300 MeV one has found that — > 1 fm is sufficient. 

Several collaborations, using dynamical quarks with pion mass down to about 300 MeV, have cal- 
culated the pion electromagnetic (EM) form factor |25], which is obtained from the matrix element 
(vr+fp/)! J^|7r+(pi)) = {pf^ + Vi ti)F-K{,'f'), where = {pj — Pi)"^ = —Q^. Recent results are shown in 
Fig. [7 |26j. Based on vector dominance, lattice data are fitted to the form F,r(Q^) = (1 + ('^^)Q^/6) ^ 
to extract the mean squared radius, which is shown in Fig. |8j As can be seen, there is an increase in 
the value of (r^) at small pion mass, m^. An accurate extraction of (r^) benefits from evaluating the 



6 



form factor at small values of accomplished by using twisted b.c. In a recent calculation, ETMC 
combined twisted b.c. and the so called 'one-end' trick to incorporate the all-to-all propagator and 
improve statistics. Using iVj = 2 simulations with twisted mass fermions at two lattice spacings and 
two volumes [27] the assessment of cut-off and volume effects was carried out. LQCD results on 
estimated in the continuum limit at pion masses in the range of 300 MeV to 500 MeV, are extrapolated 
to the physical point using NNLO chiral perturbation theory (PT). The resulting form factor is shown 
in Fig. [9] [27j and it is in agreement with experiment. 
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Figure 7: Recent LQCD results 
on as a function of [26] in 
units of the scale [2H]- 



Figure 8: The pion mean 
square radius function of 
ml with A^/ = 2 TMF. 



Figure 9: at the physi- 

cal pion mass (red band) com- 
pared to experiment (blue). 




As simulations with quark masses close to the physical value become available, the study of res- 
onances and decays of unstable particles becomes an important issue. The p-meson width has been 
studied by several groups [29]. Considering a tt+tt^ system in the / = 1-channel, the P-wave scattering 
phase shift 5ii(fc) in infinite volume is related via Liischer's relation to the energy shift in a finite box. 
Using N j = 2 TMF and considering the center of mass frame and two moving frames one extracts 



the phase shift at different values of the energy, shown in Fig. 10 From the effective range formula 



tan5ii(A;) = — / 'f f2 — Y' where k = J Eqj^/A — ml one determines mR and the coupling gp-„T, 
and then extracts the width using Tp = where = \Jm\/A — m^. The results on the width 



as a function of m^ are shown in Fig. 11 [30 
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2.3 Disconnected contributions 



Recently, progress has been made in the evaluation of the disconnected contributions to the three- 
point functions using stochastic techniques |311 132] . A case study was carried out to compare various 
stochastic methods to the exact evaluation enabled using graphics cards. For this test case, Nf = 2 
Wilson fermions simulated by the SESAM Collaboration on a volume of 16^ x 32 at m.^ ~ 750 MeV 



were used. In Figs. 12 and 13 we show results on the connected and disconnected contributions to the 
nucleon magnetic and scalar FFs, respectively. As can be seen, the disconnected contribution to the 
the magnetic FF is consistent with zero, whereas to the scalar FF is of the same order as the connected 
one. Furthermore, the scalar disconnected part converges with much fewer stochastic noise vectors as 
compared to disconnected loops contributing to the nucleon FFs, which show slow convergence [33] . A 
particularly suitable method for evaluating these fermionic loops is the truncated solver method 
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Figure 12: Connected (upper) and disconnected Figure 13: Connected (upper) and disconnected 
(lower) contributions to the nucleon magnetic FF. (lower) contributions to the nucleon scalar FF. 



3 Nucleon Generalized form factors 



High energy scattering can be formulated in terms of light-cone correlation functions. Considering 
one-particle states \p') and \p), GPDs are defined by |8l|9]: 



1 fd\ ■ 



e«'-^(p'|^(-Ar2/2)rpe 



A/2 

ig J dan-A{na) 

-A/2 



ij{Xn/2)\p) 



2 J 271 

where P = {p' + p)/2, ^ = —n ■ q/2, x is the momentum fraction, and n is a light-cone vector with 
P ■n = l. 

There are three different types of operators, depending on the choice of P . "Handbag" diagram 
Considering nucleon states these are: 



P = t/i:-^ -Un{p') 



2m 



fii/ 



N 



r = il5-^^UN{p') ii-i^H{x,i,q^) + ^^^^E{x,i,q^) 



un{p) 
un{p) 




r = n^a^'^ :^ tensor GPDs . ^ P' 

The forward proton matrix elements Fp(x,0,0), measured in deep inelastic scattering, are connected 
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to the parton distributions q{x), Aq{x), 6q{x). Expansion of the hght cone operator leads to a tower of 
local twist-2 operators ^ ^j^g forward matrix elements of which are related to moments: 



^^^1 I unpolarized 



helicity 



x")a9 = / rfxx" [Aq{x) + (-l)"Ag(x)] 



O'sg^'-^'' = i^af^^lD^' . . . iD^-^^ transversity ^ ^„ ^^^^^^ _ 

where q = q^ + q'\,A.q = qx — q^,5q = q-j + q±, and the curly brackets represent a symmetrization 
over indices and subtraction of traces. The off-diagonal matrix elements extracted from deep virtual 
Compton scattering can be written in terms of generalized form factors (GFFs), which contain both 
form factors and parton distributions: 



1 



+mod(n,2)C„+i,o(g')-gV^ ...q 

m 



un{p, s) 



and similarly for C^'^^^-'^" (in terms of Bni{q^)) and O^^'-^" (in terms of A^,, Cj„ D^J. 

In this work we will show results on the following special cases: 
• n = 1: Ordinary nucleon form factors: 

^io(g') = i^i(g') = f^dxH{x,^,q^), 5io(g') = F,{q^) = j'^dx E{x,^,q') 



Aioiq') = GA{q' 



dxH{x,^,q^), B,o{q') = Gp{q' 



dxE{x,i,q^) 



where in the case of the EM current, = ilj{x)'-f^ilj{x), the nucleon matrix element is written in the 
form uj\f{p',s') '^^Fi{q^) + ^2m^ ^^{q^) Uiy{p,s). The Dirac Fi and Pauli F2 FFs are related to the 



(2miv)' 



,F2(g2) and GmI?' 



electric and magnetic Sachs FFs via the relations: Gsiq'^) = Fi{q^) — 

Fi[q^) + F2{q^). For the axial- vector current = '?/'(x)7^75^^/'(x) the nucleon matrix element is of 

the form un{p', s') 7^75^/1(9^) + ^^Gpi^q"^) ^un{p, s), where Ga{0) = qa is the nucleon axial charge. 

• A„o(0), y4„o(0), ^^0(0) moments of parton distributions, e.g. {x)q = ^20(0) and {x)^g = ^20(0) 
are the first moments of the unpolarized and helicity distributions. Knowing these quantities one can 
evaluate the quark spin, Jg = ^[^420(0) + -820(0)] = ^AE^ + Lg and investigate the fraction of the spin 



carried by quarks and its contribution to the total spin via the nucleon spin sum rule, ^ = ^ AS„-|-Lg+ J, 



5' 



as well as the momentum fraction carried by gluons via the momentum sum rule: {x)g = 1 — ^420(0). 

In order to compare LQCD matrix elements to physical observables we need to renormalize them. 
Most collaborations carry out a non-perturbative evaluation of the vertex function r^i,(p). The renor- 
malization constants can be determined in the RI'-MOM scheme by imposing the following conditions 
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Tt[{G{p))-'G^'\p)] 



1 



Zg'Zo-TT[TMTfy\p)] 



(10) 



to extract Zg and Zo. Subtracting (9(a^)-terms perturbatively improves their determination at the 
continuum limit i34j. 
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3.1 Results on nucleon form factors 



A number of lattice QCD groups have recently produced results on nucleon FFs employing dynamical 
quark simulations with (9(a)-improved actions and lowest pion mass typically round ~ 300 MeV |12j . 
• Nucleon axial charge: The axial charge is well known experimentally. Since it can be determined 
at = there is no ambiguity associated with having to fit the Q^-dependence of a FF, such as, for 
example, in the case of the anomalous magnetic moment where one needs to fit the small Q^-dependence 
of the magnetic FF Gm- In addition, only the connected diagram shown in Fig. [5] contributes. In Fig. 14 
we show recent LQCD results using TMF, Clover fermions, domain wall fermions (DWF) and a hybrid 
action of DWF on a staggered sea, all of which are renormalized non-perturbatively. As can be seen, 
there is a nice agreement among different lattice discretizations and no significant dependence on the 
quark mass down to about = 270 MeV. 
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Figure 14: Left: LQCD results on qa using Nf = 2 TMF [35j (filled symbols, star and open square), 
Nf = 2 + 1 DWF [36] (crosses), Nf = 2 + 1 hybrid action [37j (open circles and triangle) and Nf = 2 
Clover [38] (square with cross). The experimental value is shown by the asterisk. Right: Volume 
corrected TMF results extrapolated to the continuum limit together with a chiral fit using HBxPT 
(blue band). The band bounded by the lines is the chiral fit to the TMF data shown on the left panel. 



To assess lattice artifacts and obtain a value of gA at the physical point, we use LQCD results obtained 
with TMF [35]. The volume corrected [39] data are extrapolated to a = using three lattice spacings, 
namely a = 0.089 fm, 0.070 fm and 0.056 fm. The continuum volume-corrected results are shown in 



Fig. 14 A chiral extrapolation using one-loop heavy baryon chiral perturbation theory (HBxPT) in 
the small scale expansion (SSE) [10] with three fit parameters produces a value of qa = 1.12(8) at the 
physical point, which is lower than the experimental value by about a standard deviation. The large 
error is due to the strong correlation between the A axial charge and the counter-term involved in the 
fit. In Sec. 5 we discuss the evaluation of the A axial charge within LQCD and examine the resulting 
chiral extrapolation. Fitting the uncorrected LQCD results obtained at the three /3-values produces the 
band shown by the dotted lines. This indicates that, for pion masses larger than ~ 300 MeV, volume 
and discretization errors are small compared to the uncertainty in the chiral extrapolation. 



• Nucleon form factors: Recent LQCD results on the EM isovector and axial FFs are shown in Fig. [15 
We observe a nice agreement among LQCD results, in particular for GeIQ"^) and Ga(Q^)- However, 
both GEiQ"^) and G'a(Q^) decrease with increasing Q"^ less rapidly than measured in experiment. We 
note that a good description of the Q^— dependence for both Ge^Q"^) and Gm{Q'^) is provided by a 
dipole form using the lattice-computed p— meson mass. If one uses HB^PT to one-loop, with explicit 
A degrees to perform a chiral extrapolation of the Pauli and Dirac form factors at small Q^-values 
then one qualitatively recovers the correct slope of the experimental data [41j- LQCD on Gp{Q'^) using 
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Figure 15: Left: Isovector electric and magnetic nucleon FFs at ~ 300 MeV using TMF [H] , 
DWF [42j, hybrid [37] and Clover fermions |13]. Experimental data are shown with the filled green 
circles accompanied with Kelly's parameterization shown with the dashed line. Right: Axial nucleon 
FFs. The solid line is a dipole fit to experimental data for G^IQ"^) combined with pion pole dominance 
to get the solid curve shown for GpiQ"^). 



TMF and those obtained using the hybrid action on a larger volume, differ at small where Gp{Q'^) 
increases rapidly due to the pion-pole behavior. This may indicate that volume effects are not negligible 
on FFs such as Gp{Q'^), which are strongly affected by the pion-pole. 



3.2 Results on nucleon moments 
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Figure 16: Left: Recent LQCD results on the isovector = {x)u-d- The physical point is from 

Ref. IB]. Right: Chiral fit to the continuum extrapolated TMF results using HB^PT. 

In this section we show results on the nucleon matrix element of the one-derivative operators 
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u — d'j^fj^ Du} d and m757{^ Du} u — c?757{^ Dv} d in the MS scheme at a scale /x = 2 GeV. In 
Fig. 16 we compare recent LQCD resuhs on the isovector momentum fraction using TMF [35], Clover 
fermions from QCDSF |SH], DWF from RBC-UKQCD [46J and LHPC |3Z] with different lattice spac- 
ings, and using a hybrid action by LHPC [3Z1- The moment of the helicity distribution A'^q'^ = {x)Au-Ad 
shows a similar behaviour. All collaborations, except LHPC, use non-perturbatively computed renor- 
malization constants. The m^-dependence of these moments in HB^PT [18] are given by: 
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Using = 1 GeV^ and the TMF results we obtain the band shown in Fig. 16, which yields a value 
higher than experiment. This is also true for {x)i^u-Ad- The very recent result by QCDSF [38] at 
~ 180 MeV remains higher than experiment and highlights the need to understand such deviations. 



3.3 Study of excited state contributions 

As we have demonstrated in the previous sections, there are discrepancies between LQCD results and 
well measured quantities such as qa- In this section we examine whether excited state contributions 
could be a possible origin of such discrepancies. This is done by performing a high-statistics analysis 
using Nf = 2 + 1 + 1 TMF configurations with m^r ~ 380 MeV at a = 0.08 fm. In order to probe 
effects of excited states the sink-source time separation was increased from about 1 fm used for the 
TMF results shown in the previous sections to about 1.5 fm for qa and to 1.9 fm for {x)u-d [49j. The 
insertion time of the operator from the source was fixed at 0.7 fm for qa and at 0.85 fm for {x)u-d- 
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Figure 17: LQCD results on qa (left) and {x)u-d (right) as a function of the sink-source time separation. 
The larger grey band is the the value obtained by using a sink-source separation of about 1 fm. 



From the results shown in Fig. [T7]we conclude that, for = 380 MeV, qa shows no contamination from 
excited states, whereas {x)u-d decreases by about 10% as compared to the value found for source-sink 
time separation of about 1 fm showing that excited state contributions are responsible for at least part 
of the discrepancy between LQCD and the experimental value. 



3.4 Nucleon spin 

Assuming that disconnected contributions remain small for pion masses down to about 300 MeV, we can 
evaluate the total spin carried by quarks in the nucleon considering only the connected contributions. 
In Fig. 18, we show results using Nf = 2 TMF for pion masses of 270 MeV < m.^ < 500 MeV [15] 
and results obtained using a hybrid action of DWF on a staggered sea [37]. As can be seen, both sets 
of data are in qualitative agreement and at the physical pion mass they both lead to a total spin of 
J" ~ 1/4 and J'' ~ for the u- and d-quarks, respectively. 
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Figure 18: Total spin for u- and d- quarks, for TMF (left) |l5] and for a hybrid action (right) 



4 to A transition form factors 

The nucleon is a spin-1/2 particle, and therefore its transverse charge densities do not exhibit a 
quadrupole pattern, nor do they encode any information on its shape. Such information can how- 
ever be obtained from the 7*A^ — )■ A transition charge densities. It is customary to characterize the N 
to A transition in terms of the three Jones-Scadron FFs, G\.ji, G%2 and G*(j2i denoting the magnetic 
dipole, electric quadrupole and Coulomb quadrupole transitions respectively [50] : 



Miv(P,s) (11) 



1/2 



^^''^ (^yiTr^rK^TrVi ) ' • '^^^ magnetic dipole transition FF is the dominant one and has been 



where ^ - 3 

accurately measured. Deformation is encoded in the ratios EMR and CMR, given by 
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Q+Q- G*Q2 



Ml 



AMi Gl 



(12) 



Ml 



with Q± = y (Ma ± M^Y + Q^- Dedicated experiments have yielded accurate measurements of the 
EMR and CMR excluding a zero value. This implies deformation in the N / A system as illustrated in 



Fig. 19 




Figure 19: Precise data at = o.l26 GeV^ 
strongly "suggesting" deformation in the N/A 
system [511152]. 
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Figure 20: Experimental and LQCD results on 
the magnetic dipole FF G*j^i versus Q^. 
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Figure 21: Results on EMR (left) and OMR (right). 



The to A transition FFs can be calculated within LQCD in a similar way to the FFs of the nucleon 
by evaluating the A-A^ three-point function associated with the EM current and the A and two-point 
functions and forming the ratio 



1/2 



(13) 



which yields, in the plateau region, the to A matrix element of the EM current. A number of 
techniques have been implemented to accurately extract the sub-dominant quadrupole form factors 
such as optimized sources that isolate the two quadrupoles from the dominant magnetic dipole and 
coherent sinks that increase the statistical accuracy |12]. 

Recent LQCD results are shown in Figs. 20 and 21 using a hybrid action of A'^^ = 2 -|- 1 staggered 
fermion simulations and domain wall valence quarks, as well as Nf = 2 -|- 1 DWF simulated by the 
RBC-UKQCD collaborations with lowest pion mass of about 300 MeV |53]. The slope of G%^i at low 

remains smaller than what is observed in experiment underestimating G\^^{0) i.e. one observes the 
same effect as for the nucleon form factors. Since G*^2 are underestimated at low like G\^^ 

is, taking ratios may remove some of these discrepancies. Indeed the EMR shown in Fig. 21 is in better 
agreement with experiment than the magnetic dipole FF, whereas CMR approaches the experimental 
values as the pion mass decreases. Despite the increase in statistics, the errors on the sub-dominant 
ratios when using DWF are large and to produce results at lower than 300 MeV pion masses to a 20% 
accuracy one would need to increase significantly the number of statistically independent evaluations. 

The A-A^ axial-vector matrix element (A(p', s')\A^\N{p, s)) can be written as 



' Gi{q') 



7^ + 



Gt{q'] 
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(14) 



where the dominant FFs C^(g^) and C^{q^) are analogous to the nucleon GA^q^) and Gp^q"^), respec 
tively. The A-A^ matrix element of the pseudo-scalar current is given by 



2mg{A{p',s')\P'\N{p, s))=A- 



-Uau{p',s')- Un{p,s) 



2m AT 



Using the axial Ward identity one obtains the relation: G^{q^ 

pole dominance relates to Gt^na through ^C^{q^) ~ 
Goldberger-Treiman (GT) relation 



+ ^ 

2 rn^—q^ 



1 G^NA{g^)fTTml 



(15) 



Pion 



GttNA f-K — '^TTT'nG^ 



which leads to the non-diagonal 

(16) 



In Fig. 22 we show results on the dominant axial FFs and on the ratio Gt^naI-kI'^i^nG^, which should 



be unity if the GT relation holds. This ratio approaches unity for Q > 0.5 GeV [53 
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Figure 22: The dominant axial to A FFs and the ratio of Eq. (16) using a hybrid action and DWF 



5 A form factors 

The A matrix element of the EM current (A(p', s')\j^\A{p, s)) can be written in terms of 4 FFs: 
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7^ + 



2Ma 



u A/sip, s) (17) 



with e.g. the quadrupole FF given by: Ge2 = (F* - tF;)-1{1+t) {F* - tF*), where r = -q^/{AM\). 
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Figure 23: Left: LQCD results on the A electric quadrupole FF. Contours of A (middle) and VL (right), 
with 3/2 spin projection along the x-axis. Dark colors denote small values. 

Using LQCD results on Ge2 one can obtain the transverse charge density of a A in the infinite momen- 



tum frame [HU |55]. This is shown in Fig. 23, where a A with spin 3/2 projection along the x-axis is 
elongated along the spin axis [5l|. In the same figure we also show the corresponding charge density of 
the Vl~ [56], which shows a similar deformation as the A. 

The A matrix elements of the axial-vector current {A{p', s')\A^^\A{p, s)) is given by 
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Ua{p',s' 



4Ml 



2Ma 



up{p,s) . (18) 



LQCD results on the dominant FFs gi and are shown in Fig. |24j The A axial charge is derived from 
gi{0). The A matrix element of the pseudo-scalar current (A(p', s')|P^|A(p, s)) is given by 



-Ua{p',s') 



2mg(m2 — q^) 



upiP, s), 



(19) 



in terms of two ttAA pseudo-scalar FFs [57], of which G^^aa is shown in Fig. 24 
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Figure 24: LQCD results on the dominant axial gi and g^, and pseudo-scalar GttAA A FFs in the 
quenched theory and using a hybrid action of DWF on a staggered sea 
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Having a set of lattice results for the axial nucleon charge [35], the axial A^ — A transition FF, C^(0) |53] 
and the A axial charge, allows us to perform a combined fit to all three quantities using HB^PT in the 
SSE scheme [SHlEniEn]. 

For gA we use results obtained with TMF since, as dis- 
cussed in Sec. 3, different discretization schemes are in 
agreement. The A axial charge and C^(0) were both 
evaluated using the same fermionic action, namely Nf = 
2 + 1 DWF and a hybrid action of DWF on staggered 
sea. The combined fit has seven free parameters: the 
three axial coupling constants of the nucleon, the iV to A 
and the A, the three coefficients of the m^-terms in each 
chiral expansion of gA, ^5^(0) and Gaa = —3(71(0) and 
a constant entering the chiral expression of C^(0) [5^ . 
As can be seen in Fig. 25, lattice data for these three 
observables show a weak pion mass dependence within 
the mass range considered. The resulting fits are shown 
by the bands. The value of gA, as computed in all recent 
lattice studies, is underestimated and this combined fit 
does not provide a possible resolution to this puzzle. Hav- 
ing lattice results at pion masses below 300 MeV will be 
essential to check the validity of these chiral expansions. 
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Figure 25: Combined chiral fits: (a) gA with Nf = 2 TMF [3^; (b) Real part of axial N to A transition 
FF (75(0) [53|; (c) Real part of A axial charge Gaa = — 3(y'i(0) ET 
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7 Conclusions 

We have shown that lattice QCD successfully reproduces the low-lying baryon spectrum using different 
discretization schemes. Understanding nucleon structure within LQCD is a fundamental issue that is 
being addressed by a number of lattice collaborations [121 EH HSl [Ml EH US] . Using similar techniques 
one can study transitions and resonant properties and in this work we have reviewed results on the 
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p-meson width, on the to A transition FFs and the A FFs. The latter are difficuh to measure 
experimentally and therefore LQCD provides valuable input on these quantities. The study of the 
complete N / A system enables one to extract the axial couplings from a combined chiral fit to LQCD 
results on the nucleon and A axial charges and the axial N to A form factor C^{Q) [57j. Performing 
such a fit to lattice results in the pion mass range from 500 MeV to 300 MeV, does not reproduce the 
experimental value of qa- Gauge configurations with pion masses below 200 MeV are becoming available 
enabling the evaluation of these key observables at near physical quark mass parameters. Using these 
simulations, combined with a detailed study of lattice systematics, are expected to shed light on the 
origin of the observed discrepancies. 
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